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ABSTRACT 

We present a sample of 824 solar and late-type stars with X-ray luminosities and rotation periods. 
This is used to study the relationship between rotation and stellar activity and derive a new estimate 
of the convective turnover time. From an unbiased subset of this sample the power law slope of the 
unsaturated regime, Lx/ Lboi oc Ro^ , is fit as = — 2.70±0.13. This is inconsistent with the canonical 
P = —2 slope to a confidence of 5a, and argues for an additional term in the dynamo number equation. 
From a simple scaling analysis this implies Aft/ft oc f2 - 7 , i.e. the differential rotation of solar-type 
stars gradually declines as they spin down. Super-saturation is observed for the fastest rotators in our 
sample and its parametric dependencies are explored. Significant correlations are found with both 
the corotation radius and the excess polar updraft, the latter theory providing a stronger dependence 
and being supported by other observations. We estimate mass-dependent empirical thresholds for 
saturation and super-saturation and map out three regimes of coronal emission. Late F-type stars 
are shown never to pass through the saturated regime, passing straight from super-saturated to un- 
saturated X-ray emission. The theoretical threshold for coronal stripping is shown to be significantly 
different from the empirical saturation threshold (Ro < 0.13), suggesting it is not responsible. Instead 
we suggest that a different dynamo configuration is at work in stars with saturated coronal emission. 
This is supported by a correlation between the empirical saturation threshold and the time when stars 
transition between convective and interface sequences in rotational spin-down models. 
Subject headings: stars: activity - X-rays: stars - stars: evolution - stars: late-type - stars: coronae - 
stars: magnetic fields - stars: rotation 
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1. INTRODUCTION 

Stars across the Hertzsprung-Russell diagram are 
known to emit X-rays with only a few exceptions. 
The hottest and most massive stars emit X-rays that 
arise from either s mall-scale shocks in their winds 
(jLucv fc White 1980h or collisi ons between the wind 
and circumstellar material (e.g. lParkin et al. 2009D . So- 
lar and late-type stars emit X-rays from a magneti- 
cally confined plasma at temperat ures of several mi l- 
lion Kelvin known as a corona (jVaiana et al. 19811 ). 
This high-temperature rarefied thermal plasma, first 
observed on the Sun, has been detected from nearly 
all classes of low-mass stars at leve ls of Lx/Lboi ~ 
1CT 8 - 10" 3 (|Schmitt fc Liefke 20041) . The most X- 
ray luminous of single, low-mass stars have been found 
to be the youngest, where X-ra y luminosity levels o f 
Lx/Lboi ~ 10 are common (e.g. iTelleschi et al. 2007D . 
though close or interacting binary stars such as RS CVns 
may also reach such levels. Older stars have typ- 
ically lower X-ray luminosities, reaching down to 
LxlL bo i ~ 10~ 8 - 10~ 4 (e.g. iFeigelson et al. 2004 
IWright et al. 20 10). This significant decrease in stel- 
lar X-ray luminosity of several orders of magnitude be- 
tween the zero age main sequence and the solar age is 
widely used to identify young sta rs efficiently against 
the Galactic field pop ulation (e.g. IGetman et al. 2005 
IWright fc Drake 2009f ). 

Coronal X-ray emission, as well as other indica- 
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tors of stellar a ctivity (e.g. Ca II or Ha emission, 
IWest et al. 20081) . is believed to be driven by the stel- 
lar magnetic dynamo, which itself is thought to be 
driven by differential rotat i on in the ste llar interior (e.g. 
iParker 19551 IWilson 196(1 : Kraft 19671) . a phenomena 
that has be en confirmed in th e Sun through helioseis- 
mology (e.g. iDuvall et al. 19841 ). The observed decrease 
in X-ray emission over the lifetime of the star can there- 
fore be attributed to the rotational spin-down of the star, 
which is driv en by mass-loss t hrough a magnetized stel- 
lar wind (e.g. lSkumanich 19721 ) . Because the magnetized 
stellar wind i s also believed to be d riven by the stellar 
dynamo fe.g. ICranmer fc Saar 20flT) . rotation and mag- 
netic activity effectively operate in a feedback loop with 
rotation, stellar activity, and the stellar wind decreas- 
ing with age. Since angular momentum loss controls the 
evolution of close binaries that are believed to be the 
progenitors of cataclysmic variables, novae, and Type la 
supcrnovae, understanding these dependencies has far- 
reaching implications. 

The relationship between stellar rotation and tracers of 
magnetic activity is an important probe of the physical 
dynamo process, with chromospheric and coronal trac- 
ers increasing monotonically with increasing rotational 
velocity for periods exceeding a few days. A relation- 
shi p between rotation and activity was first quantified 
by iPallavicini et al. fl98il) who found that X-ray lumi- 
nosity scaled as Lx oc (usini) 1,9 , providing the first ev- 
idence for the dynamo-induced nature of stellar coro- 
nal activity. For very fast rotators the relationship was 
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found to break down ijMicela et al. 19851) , with X-ray lu- 
minosity reaching a saturation level of approximatel y 

log Lx/L bol 3 dVilhu 1984t iVilhu fc Walter 19871) . 

independent of spectral-type. This saturation level is 
reached at a rotation period that increases toward later 
spectral t ypes (increasing with decreasing bolometric lu- 
minosity, IPizzolato et al. 20031) . It is unclear whether 
this satur ation is caus ed by a saturation of the dynamo 
itself (e.g. IVilhu 1984D . a saturation o f the filling factor 
of active regions on the stellar surface (|Vilhu 1 984) , or a 
centrifugal s tripping of the corona ca used by the high ro- 
tation rates (jJardine fc Unruh 19991 ). However, once sat- 
uration occurs the X-ray emissi on becomes a functio n of 
only the bolometric luminosity (|Pizzolato et al. 20031 ) , or 
effectively the mass, color or radius of the main-sequence 
star. 

In the non-saturated regime, the two influences on the 
efficiency of the m agnetic dynamo were combined by 
iNoyes et al. (19841 ) into a single parameter, the Rossby 
number, Rq = P ro t/T, the ratio of the stellar rota- 
tion period, P ro t, and the mass-dependent convective 
turnover time, r. This quantity has proven to be 
an effective parameter of the stellar magnetic dynamo, 
increasing towar d lower masses with the efficiency of 
the dynamo ( e.g. IMicela et al. 1981 iMaggio et al. 19871; 
IStepien 19941 : iRandich 20001 ). Despite this work there is 
yet to be a satisfactory dynamo theory that can explain 
both the solar dynamo and that of rapidly rotating stars 
(e.g. [Weiss 200i| iBrandenburg et al. 20111 ) and the con- 
tinued lack of a sufficiently large and unbiased sample 
has no doubt contributed to this. 

The paucity of stellar samples with which to study 
the rotation-activity relationship has mainly been due 
to the difficulty of measuring accurate stellar rotation 
periods, which require multiple deep observations over 
long baselines. This has led to the use of projected rota- 
tional velocities as a substitute, which are influenced by 
the uncertainties of estimated stellar radii and unknown 
inclination angles. The recent increase in measured ro- 
tation periods (e.g. lHartman et al. 20 101 increased the 
number of Pleiades stars with measured periods by a fac- 
tor of five) for many thousands of stars in open clusters 
of known age is overcoming this problem and it is likely 
that we will henceforth be limited by the availability of 
deep X-ray observations for such stars. 

In this work, we combine new measurements of photo- 
metric rotation periods for a large number of field and 
cluster stars with archival X-ray observations to produce 
the largest existing sample of stars with photometric ro- 
tation periods and X-ray luminosities (Section 2). This 
sample is then used in Section 3 to study and charac- 
terize the rotation - activity relationship in detail and 
to probe the stellar magnetic dynamo responsible for it. 
This allows us in Section 4 to trace out the X-ray evo- 
lution of low-mass stars as a function of rotation period, 
which is a good proxy for age. Finally, in Section 5 this 
sample is used to derive a new empirical measure of the 
mass-dependent convective turnover time. 

2. COMPILATION OF THE SAMPLE 

To study the relationship between rotation and activity 
a sample was compiled from the literature by searching 
for stars with measurements of both rotation periods and 
X-ray luminosities. Only photometrically-determined ro- 



tation periods were included, discarding all rotation ve- 
locity measurements and upper limits, and only stars 
with significant X-ray detections were used, discarding 
all sources with only upper limits. This choice reduces 
the sample size available and also has the potential to in- 
troduce an X-ray luminosity bias in our results. However 
this greatly simplifies the following analysis, particularly 
in the light of the large variety of sources used to com- 
pile this sample, the different techniques used to calcu- 
late upper limits by different authors, and the potential 
incompletenesses in upper limits present in each sample. 
The inherent biases that will exist in this sample will be 
discussed and addressed later. 

2.1. Rotation periods and X-ray fluxes 

The recent study o f the activity-rotation relation by 
IPizzolato et al. (20031 ) provided the starting point for the 
catalog. From their work 102 cluster stars (excluding 
Pleiades members, which were compiled separately) and 
47 field stars were used, excluding all sources with up- 
per limits, as well as a number of stars whose rotation 
periods were inferred indirectly from chromospheric ac- 
tivity levels (with the exception of a Ccntauri B, for 
which we use the rota tion period and mea n X-ray lu- 
minosity presented by iDeWarf et al. (2010t )). The ma- 
jority of these sources are G and K stars with the M- 
type stars confined to the field star sample because of 
the lack of rotational periods available for low mass stars 
in clusters at the time. A further 28 Hyades members 
were intr oduced by cross-matc hing recent rotation peri- 
ods from iDelorme et al. (201 1[) with X-ray lum inosities 
from IStern et al. (19941) and lStern et al. (19951 ) . 

The first expansion of the sample was based on the 
rec ent measurement of photometric rotation periods 
by lHartman et al. (20101 ) for Pleiades membe rs based 
on the membership list of IStauffer et al. (20071 ). These 
were cross-matched with X-ray flux measurements of 
Pleiades sta rs from the Einstei n and ROSAT observa - 
tions from IMicela et al. (1990D. IStauffer et al (19941) . 
IMicela et al. (1996D . and IMicela et al. (19991) . an d 
XMM-Newton observations from iBriggs fc Pve ( 2003). 
again discarding all upper limits. Where multiple 
measurements of the X-ray flux exist for a single source, 
that with the lowest fractional uncertainty was used. 
Cross-matching 391 X-ray sources with the 383 sources 
with photometric periods lead to a sample of 146 
Pleiades members with both X-ray luminosities and 
rotation periods, the largest such sample for a single 
cluster. 

This sample was then complemented with 83 stars 
from the open clusters NGC 2516 an d NGC 2547 using , 
respe ctively, rotation p eriods from llrwin et al. (2007 ) 
and | Irwin et al. (20081) . and X -ray flux es from 
iPillitteri et al. (20061 ) and Ueffries et al. (201l1) . An ad- 
ditional 20 stars were added from the open cluster Prae - 
sepe using rotation per iods from IDelorme et al. (201 ID 
and IScholz et al. (201 ll) and X -ray luminosities from 
IRandich fc Schmitt (19951 ) and iFranciosini et al. (2001 . 
The sample was then further extended using data 
from rotation period surveys of field stars combined 
with X-ray fluxes from the R OSAT All- Sky Survey . 
This included 218 stars from lHartman et al. (201ll ) , 
23 stars from iKiraga fc Stepien (20071). 6 stars from 
IXing et al. (2007| ). 8 stars from iBouvier et al. (19971 ) 
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and 79 stars were added fr om the compilation of 
lMamajek~ fc Hillen brand (20081 ). This sample also in- 
cludes the Sun usi ng the values of log Lx / L^oi = —6.24 
([Judge et al. 20031) 1 and P rot = 26.09 days 
( Donahue et al. 19961 ) . Finally, rotation periods for 
65 stars were taken from observations as part of the 
FEPS (Formation and Evolution of Planetary Systems, 
see Appendix A) program. 

Known classical T-Tauri stars with Ha emission, sig- 
naling the presence of accretion and therefore a cir- 
cumstellar disk, were excluded because of the complica- 
tions induced by X-ray emission from accretion and disk- 
locking on the rotation period. Pre-MS stars (specifically 
those with ages <10 Myrs) were also excluded because 
of potential differences in their internal structure as a 
function of either mass or effective temperature. Us- 
ing the catalog of X-ray vari able sources presented by 
iFuhrmeister fc Schmitt ?2003l) . we removed all ROSAT 
sources that had been observed to flare during the ob- 
servation to lessen the influence of short-duration flares 
on the resulting X-ray luminosities. In total the result- 
ing sample includes 824 stars, all with photometrically- 
measured rotation periods and well-constrained X-ray 
flux measurements. 

2.2. Homogenizing the sample 

The majority of previous work on the rotation-activity 
relation has used the B — V color as the proxy for effec- 
tive temperature. This was due primarily to the ubiq- 
uity of suitable photometry for the FGK stars that made 
up these studies and the lack of photometry in other 
bands. However, the increased interest in M-type stars, 
for which the B — V color is a poor temperature diag- 
nostic, and the availability of near-IR photometry sug- 
gests an alternative color as an appropriate proxy. The 
V — K s color, derived from now readily available optical 
and near-IR photometry offers a sufficiently large base- 
line such that uncertainties in either individual magni- 
tude do not greatly affect the derived color. Comple- 
mentar y photometry was ther efore sought from the Hip- 
parcos (|Perrvman e t al. 199 3) and Two Micron All Sky 
Survey (2MASS, iCutri et al. 20031 ) catalogs and used to 
derive V — K s colors for all stars from the available pho- 
tometric colors with the widest baseline, using the empir- 
ical tabulations of V — K s versus other co lors and spec- 
tral type presented bv lPecaut et al. (20111 ). 797 sources 
have V and near-IR photometry (2MASS J, H, or K s 
band observations with the lowest photometric uncer- 
tainty and photometric quality flags of A, B, or C), 17 
sources had V — I colors from which V — K s colors were 
derived, and 10 sources had only B — V photometry, 
due primarily to their proximity and the saturation of 
near-IR photometry. The literature was then searched 
for suitable spectral types for these sources (see Ap- 
pendix B for references) and these were used to derive 
the intrinsic V — K s color (using the empirical tables of 
iPecaut et al. 20ilT) for sources where the spectra would 
provide a more accurate diagnosis of the underlying pho- 
tospheric temperature than photometry could 2 . In total 

1 This value is larger than that of Peres et al. (2000) by a factor 
two, a factor that would only exacerbates the poor fit of the solar 
data point compared to the fit in Section |3. II 

2 Given a typical classification uncertainty of ±1 subtype we 
estimate that spectra provide a more accurate classification than 



TABLE 1 

List of cluster and field star samples used in this work. 



Cluster / nela 
star sample 


Age 
(Myrsj 


Distance 

(P c ) 


Extinction 

(A \ 

\ A V) 


Number 
of stars 






A (17 


U. 10 


u)t 


IC 2602 6 


46 


149 


0.09 


28 


IC 2391 


50 


145 


0.03 


13 


a Persei c 


85 


172 


0.27 


10 


Pleiades^ 


125 


133 


0.12 


146 


NGC 2516= 


150 


376 


0.36 


14 


Praesepe 


600 


182 


0.08 


20 


Hyades 


700 


46.5 


0.0 


49 


Bouvier 1997 




143 - 264 




8 


Pizzolato 2003 




1.3 - 172 




47 


Kiraga 2007 




3.0 - 50 




23 


Xing 2007 




44 - 207 




6 


Mamajek 2008 




4 - 102 




79 


FEPS program 




19 - 102 




64 


Hartman 2011 




1.5 - 243 




218 


Total 








824 



Cluster distances and extinctions arc taken from the Hipparcos 

parallax measurements dvan Le cuwcn 2009, and references 
therein), with the exception of the Pl eiades, for which a distance 
of 133 pc is used (Stauffcr et al. 1998), and the relatively distant 
clusters NGC 2516 and NGC 2547, which are poorly sampled by 
Hipparcos. Clust er ages are cited individ ually. 

Reference s. ( a)IMavne fc Navlor (20011; (b) 

IDobbie et al. (2010); (c)IBarrado y Navascues et al. (20041 ); (d) 
IStauffer et al. (19981 ); (e) I Jeffries et al. (20011 ). 



225 sources had spectral types in the literature that were 
used to derive dereddened V — K s colors, including all 
but 4 of the stars that lacked V — K s photometry. This 
also included all field stars more distant than 75 pc for 
which interstellar reddening might have been a concern 
for purely photometric data. The colors of cluster stars 
were dereddened using the best available measurements 
of cluster extinction (Table [J) to provide the most accu- 
rate available effective temperature proxy for each star 
in the sample. 

Using the V — K s color as a proxy for effective tempera- 
ture, we derived bolometr ic corrections (BC) for all sta rs 
using the tabulation of iKenvon fc Hartmann f! 995) 3 . 
Distanc es for field stars were calculated usin g Hip- 
parcos ([Perrvman et al. 1997t Ivan Leeuwen 2007t ) paral- 
laxes, considering only those with good quality measure- 
ments, i.e. relative uncertainties < 20%. For the 559 
stars with known distances (378 cluster and 181 field 
stars) we combined these distances with BCs and U-band 
magnitudes to derive absolute b olometric lumino sities. 
Combining the isochrones from ISiess et al. (2000T ) and 
the most recent estimates for the ages of each cluster 
(Table [T]), or assuming an age of 1 Gyr for field stars, we 
used the bolometric luminosities to derive stellar masses, 
radii, and effective temperatures for each star. 

For the field stars for which distances were not known 
the V — K s color was used to derive the effective 
temperatur e using the calibrated eff ective temperature 
scale from ICasagrande et al. (20 08. for M stars) and 
ICasagrande et al. (20101 for FGK stars), with a lin- 
ear interpolation between the two tabulations. The 

photometry for all F and G-type stars, for stars known to be in 
binary systems, and for any source where the B — V color was 
used. For late- type stars where either the V — K B or V — I color was 
available we consider the photometry to have a smaller uncertainty 
than a typical spectral classification. 

3 Conversions provided by Bcsscll (2005) were used to convert 
the V — K color of this work onto the 2MASS-based V — K s color. 
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effective tempera ture was then fitted to the 1 Gyr 
iSiess et al. (2 000) main sequence isochrone to derive stel- 
lar masses, radii, bolometric luminosities, and therefore 
distances for all these stars. A comparison of the differ- 
ent color-T e // tabulations suggests that the systematic 
uncertainty in the derived stellar mass using this method 
is potentially as much as 10% for FGK stars, and 25% 
for M-type stars, though we are confident we have cho- 
sen the most recent and reliable color-T e / f conversion for 
this use. For comparison, the derived mass uncertainty 
based on the color-BC conversion is < 2% for all stars 
with V — K a < 6 (the majority of stars in our sample), 
excluding uncertainties in the distance used, the age of 
the clust er or systematic errors in the evolutionary tracks 
(see e.g. iHillenbrand fc White 200l for a discussion of 
the accuracy of these tracks). 

Using the derived distances, X-ray luminosities and X- 
ray-to-bolometric luminosity ratios were determined for 
all sources. For the majority of stars X-ray fluxes were 
available in the literature (based on X-ray spectral fits, 
count rates or hardness ratios), though for a small num- 
ber of stars X-ray fluxes were either not available or had 
been calculated using a constant count-rate to flux con- 
version (i.e. one that does not take into account the 
temperature of the underlying spectrum through hard- 
ness ratios or spectral fits). This was mainly the case for 
some of the early ROSAT papers such as those study- 
ing the Hyades. For these sources w e adopted the con- 
version between ROSAT count ra tes rtVoges et al. 19991 ) 
and hardness ratios published bv lFleming et al. (1995T ). 
For comparison with the majority of literature X-ray 
observations of open clusters we have converted all 
X-ray luminosities to the ROSAT 0.1-2.4 keV energy 
band. This wa s done using PIMMS 4 , as suming a ther- 
mal spectrum ()Ravmond fc Smith 19 77) and a hydro- 
gen column density converted from the visual extinction 
(jRvter 1996( 1. The temperature of the thermal spec- 
trum for the conversion was chosen based on the ac- 
tivity level of the star, Lx/Li, i, and the observed cor- 
relation between activity leyel and plasma temp erature 
(see e.g- lTerieschi et al. 2005HGudel fc Naze 2009h 5 . Un- 
certainties in these luminosities due to uncertainties in 
the hydrogen column density and plasma temperature 
are ~5% and ~9% for dispersions of A log Nh = 0.5 
and AkT = 0.5 k eV, as appropriate for the Pleiades 
(|Gagne et al. 19 95). For reference with other works the 
conversion factors to go from the ROSAT band to the 
Chandra 0.5-8.0 keV or XMM-Newton 0.3-4.5 keV bands 
are 0.676 and 0.797, respectively, for an active star with 
kT ~ 1 keV. 

2.3. The final catalog 

The compilation of this data led to a catalog of 
824 stars with rotation periods and X-ray luminosities 
that represents a significant increase over those used 
for previous stu dies of the act ivity-rotation relation 
bv lStepien (1994L ~70 stars) and iPizzolato et al. (200l 
^250 stars). Table [1] provides the details of all of these 

4 The Chandra Portable Interactive Multi-Mission Simulator, 
http: / /cxc. harvard.edu/toolkit/pimms.jsp 

We calibrated thi s relation using the investigation of solar 
analogs by Tellcs chi et al. (20051 ) , but convert the dependency from 
the X-ray luminosity to the X-ray luminosity ratio as a more ap- 
propriate diagnostic of the activity level of the star. 
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Fig. 1. — Cumulative distribution of stars in our sample as a 
function of their V — K color (top), used in this work as the ob- 
servable proxy for effective t emperature, and s tellar mass (bottom) , 
derived from the models of Sicss et al. (2000). Reference spectral 
types are indicated. 

samples and the values used to derive the stellar param- 
eters. Figure [T] shows the distribution of the sample as 
a function of V — K s and mass. In V — K s color space 
the sample is approximately equally distributed across 
the range 1.5 - 5.0 (G2 - M4), with ~30 stars per spec- 
tral subclass, dropping to ~10 stars per subclass beyond 
these limits. Uncertainties of 10% in the source or cluster 
distances translate to uncertainties <10% in the resulting 
masses and uncertainties of ~20% in the resulting X-ray 
luminosities. However, distance uncertainties in the X- 
ray to bolometric luminosity ratios cancel out, and only 
uncertainties in the stellar mass become relevant. These 
uncertainties are generally small compared to uncertain- 
ties in the measured X-ray luminosities, which are found 
to have a typical uncertainty of ~0.2 dex. 

3. X-RAY EMISSION VS. STELLAR ROTATION 

In this section the accumulated data is used to study 
the relation between the level of X-ray activity, Lx , and 
the stellar rotation period, P, both of which vary by 
many orders of magnitude across the sample. These two 
parameters ar e thought to be connected by the magnetic 
dynamo (e.g. P arker 19551) . but the form of this rela- 
tion and influence of spectral type is poorly understood. 
Because of this it has become usual to scale both ob- 
servable parameters by function s of the stellar mass that 
allow a more useful comparison. iPallavicini et al. (1981f ) 
represented the level of activity with the X-ray to bolo- 
metric luminosity ratio, Rx = Lx/Lboi- Following 
iNoves et al. (1984 ) the rotation rate is represented with 
the Rossby number, Ro — P/t, the ratio of the stellar 
rotation period to the convective turnover time, r. The 
objective of these transformations is to convert the ob- 
servable quantities into those that represent the parame- 
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ters and products of the stellar dynamo. Figure [2] shows 
the X-ray luminosity ratio as a function of both the ro- 
tation period and the Rossby number and it is clear that 
the latter parameterization greatly reduces the scatter in 
the unsaturated regime. 

The color-dependent convective turnover time, 
r, cannot be measured directly but can either 
be derived from mo dels of stellar interiors (e.g. 
IKim fc Demargue 1 996) or be empirically estimated 
(e.g. iNoves et al. 1984ft. The semi-em pirical deter- 
mination provided by INoves et al. (19841) is the most 
well-used in the literature, but was derived from a small 
sa mple that only extends to B — V = 1.4, and as noted 
by iPizzolato et al. (20031 ) is based on only 5 points 
redder tha n B — V = 1.0. T he empirically-derived 
values from IPizzolato et al. (20031 ) offer an improvement 
on this, being based on a sample that extends to 
B — V = 2.0, with a greater coverage of low-mass stars, 
and it is therefore chosen for this work. We use an 
empirically-de rived conversion between V — K s and 
B-V colors (jPecaut et al. 20111 ) to estimate convective 
turnover times for each star. We note that the convective 
turn over time is known to v ary over the life of a star 
(e.g. IKim fc Demargue 19961 ) and may vary significantly 
in the pre-MS phase, potentially causing a systematic 
error in the analysis of the rotation-activity relation 
for such stars. In Sec tion [5] we follow the method 
of IPizzolato et al. (20031 ) to derive a new empirical 
estimate of r based on the larger sample used in this 
work. 

Figure [5] shows Rx as a function of Ro for all the stars 
in our sample, clearly demonstrating that the X-ray lu- 
minosity ratio increases with decreasing rotation period 
(or increasing rotation rate), as expected from qualita- 
tive arguments based on the a — cj-type shell dynamo 
theory. As noted by many previous observers, the X-ray 
emission level appears to saturate at the highest rota- 
tion rates, reaching an approximately constant level of 
Rx ~ 10 -3 . This effect clearly divides the rotation- 
activity relation into two regimes: a saturated regime at 
high rotation rates, and an unsaturated regime for slow 
rotators. The transition between these two regimes can 
be seen to occur at Ro ^0.1. 

In the discussion that follows, the rotation-activity re- 
lation is divided into these two regimes in an attempt to 
reveal the different physical processes at work. Approx- 
imately linear relations in log-log space are immediately 
apparent from this figure, despite a spread in both Rx 
and Ro. This spread is likely to be caused by a number 
of different factors arising from the necessarily simplified 
analysis of the data performed here. Both the X-ray lu- 
minosity and the photometric rotation period are approx- 
imations of the real values due to the methods employed 
in measuring them. X-ray luminosities will vary over 
the course of a stellar activity cycle, by up to an orde r 
of magnitude in the case of our Sun ()Peres et al. 2 000). 
and will vary on shorter timescales due to the influence 
of coronal flares. Furthermore, the measured rotation 
period will be a latitudinal mean due to the unresolved 
nature of the stellar disk, and may also vary over the 
course of an activity cycle due to the varying starspot 
coverage. 

A further source of uncertainty is due to the effects of 
unresolved binaries in this sample. Photometrically, the 



presence of a lower mass binary companion will cause 
the star to appear redder and more luminous, result- 
ing in a lower mass estimate and a smaller photomet- 
ric distance. Distance errors cancel out when deriving 
the X-ray-to-bolometric luminosity ratio, but uncertain- 
ties in the mass, and therefore spectral type, can affect 
a number of the parameters and quantities derived for 
each star. Literature spectral types were obtained for 
225 of the 824 stars in the sample including 163 of the 
168 binaries in our sample, reducing the influence of this 
effect. Another influence of binary companions is the 
effect on the X-ray emission, the presence of a close sec- 
ondary potentially elevating the level of emission due to 
tidal coupling or their orbits. The 168 binaries in our 
samples are highlighted in Figure [2] and are concentrated 
in the unsaturated regime (which is dominated by older 
field stars). Their X-ray levels, in both the saturated and 
unsaturated regimes, all lie within the la spreads on the 
fits determined here, with no trend for cither elevated 
or reduced X-ray emission. We therefore include these 
sources in our subsequent study of the rotation-activity 
relation without fear that their presence will bias the re- 
sults. 

3.1. Dynamo efficiency in the unsaturated regime 

The unsaturated regime in the rotation-activity re- 
lation is believed to probe the efficiency of the stellar 
dynamo in heating the corona. iPallavicini et al. (198 J ) 
found that the X-ray luminosity of solar- and late-type 
stars scales with projected rotational velocity, to the first 
order, Lx oc (wsini) 2 . This relationship has since been 
investigated by many authors, more recently in the form 
of the Rx - Ro relationship, using photometric rotation 
periods as well as making the distinction between stars 
with saturated and unsaturated X-ray emission. The low 
level of scatter in the Rx — R o diagram has been in- 
terpreted by many a uthors (e.g. iMontesinos et al. 20011 : 
IPizzolato et al. 20031 ) as evidence for a strong underlying 
physical relationship. 

To parameterize this relationship the sample in Fig- 
ure [2] was fitted with a two-part function of the form 

R _ J Rxsat A Ro < Ro sat , v 

Ux - \CRoP if Ro > Ro sat [L> 

where the parameters Rxsat, Ro sa t, and ft (the X-ray 
luminosity ratio and Rossby number at saturation, and 
the power-law index) were varied to find the best fit us- 
ing a ^-minimization technique. Uncertainties on these 
quantities were then determined using a bootstrapping 
approach, iterating 1000 times and finding the standard 
deviation of each parameter. The parameters Rxsat and 
Ro S at will be discussed in more detail in Section 13.21 
and were found to be log Rxsat = —3.13 ± 0.08 and 
Rosat = 0.13 ± 0.02 for the best-fitting model, which 
has an rms scatter of ^0.3 dex in Rx- The best- 
fitting slope to the unsaturated regime was found to be 
(3 = —2.18 ± 0.16, slightly steeper than the canonical 
value of /? ~ —2. This fit, as shown in Figure [21 over- 
predicts the Sun's mean X-ray luminosity by a factor of 
2-3. 

An alternative approach is to fit the slope of the un- 
saturated regime and the saturation level separately. 
We fitted a simple power law of the form logRx — 
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Fig. 2. — X-ray to bolometric luminosity ratio plotted against rotation period (left panel) and the Rossby number, Ro = P ro t/T (right 
panel), for all stars in our sample with X-ray luminosities and photometric rotation periods. Stars known to be binaries are shown as plus 
symbols, and the Sun is indicated with a solar symbol. The best-fitting saturated and non-saturated activity— rotation relations described 
in the text are shown as a dashed red line in the right-hand panel. 



log C+/3 log Ro to all stars with Ro > 0.2 using the differ - 
ent types of linear regression fits in llsobe et al. (1990l ). 
We find a good agreement between the slopes derived 
from these different fits, suggesting that the fits are all 
fairly linear in the log Rx - log Ro plane. We favor 
the Ordinary Least Squares (OLS) bisector since the 
objective of the fit is to estimate the underlying func- 
tio nal relation betwe en the variables, as recommended 
bv llsobe et al. (1990| ). and this method also factors in the 
scatter of the line in both variables. The fit gives a slope 
of j3 = -2.55 ± 0.15 (valid in the range 0.2 < Ro < 3, 
or —3.75 > \og(Lx /Lboi) > —6.3), significantly steeper 
than both the canonical value and that foun d from our 
two-part fit. Mamajck & Hillenbrand (2008) fit a log- 
linear function to the Rx — Ro distribution, with the 
goal of empirically deriving a correlation that would al- 
low age estimates to be derived from X-ray luminosities 
(via rotation periods). Their fit does not connect with 
the level of saturated X-ray emission for very fast rota- 
tors but, as they note, it offers a good fit to many of the 
slow rotators such as the Sun. 

The sample used here suffers from a number of biases 
due to the selection of only sources with measured X- 
ray fluxes and photometric rotation periods. While the 
biases stemming from the detectability of rotation pe- 
riods are myriad and complex, the luminosity bias in- 
duced by only using sources with measured X-ray fluxes 
is clear. This bias will be most prominent in the unsat- 
urated regime where X-ray luminosity ratios may reach 
as low as ~ 10~ 7 or lower. This sample could there- 
fore be missing some of the faintest sources at a given 
Rossby number, possibly resulting in a larger spread in 
the Rx — Ro diagram than is currently observed. Such 
a spread could easily be induced by the increased am- 
plitude of stellar coronal cy cles that has been sugg ested 
to occur as stars age (e.g. iMicela fc Marino 20031 ) . At 
the largest Rossby numbers it is likely that many of the 
faintest X-ray sources are not detected, inducing a strong 
bias in our sample that will affect the fits derived here. 



3.1.1. Probing the dynamo efficiency with an X-ray 
unbiased sample 

To overcome the biases in our large sample we have 
attempted to compile from within our sample a smaller, 
X-ray unbiased sample that covers a large range in X- 
ray luminosity ratios and rotation periods. For this we 
use the list of 36 M t. Wilson stars with ro tation periods 
from the study by iDonahue et al. (T9 96). all of which 
were detected by ROSAT and therefore do not suffer 
from X-ray luminosity biases. These 36 stars are the 
subsamplc of their entire sample of 100 observed stars 
with measurable rotation periods over five or more sea- 
sons. The authors discuss a number of possible biases in 
their sample resulting from effects such as active region 
growth and decay, multiple active regions, and latitudi- 
nal bands. They conclude that the resulting biases affect 
only AP, not the period itself, and are either small or act 
to reduce the measured value of AP. Therefore we be- 
lieve that this sample of 36 stars with measured rotation 
periods and and X-ray luminosities is free from the ma- 
jority of biases. These stars were included in our sa mple 
as part of the compilations of | Pizzolato et al. (20031) and 
iMamaiek fc Hillenbrand (20081 ). and in Figure[3]we show 
their distribution in the Rx~Ro diagram, all of which fall 
in the unsaturated regime of coronal emission. 

We fitted a simple single-part power law of the form 
logPx = log C + log Ro to these points, using an OLS 
bisector (jlsobe et al. 19901 ) . though the slopes derived 
from all the different fitting methods are in good agree- 
ment. The fit gives a slope of/3 = —2. 70 ±0.13 (valid in 
the range 0.3 < Ro < 3, or —4 > \og(Lx /L^oi) > —6.3), 
steeper than that found from our larger sample, in agree- 
ment with our predictions of the uncertainties induced by 
the biases of that sample. This slope is even steeper than 
the canonical value of = —2 as well as the slope found 
by iPallavicini et al. (19811) of = -1.9 ± 0.5, though 
their use of projected rotational velocities instead of ro- 
tation periods represents a different relationship than 
that fitte d here. However th is slope is in good agree- 
ment with lGiidel et al. (19971 ) who derive a similar slope, 
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Fig. 3. — X-ray to bolometric luminosity ratio versus Rossby 
number for our unbiased sample of 36 stars with unsaturated X- 
ray emission. The Sun is indicated with a solar symbol. The log-log 
OLS bisector fit, = —2.70, is shown as a dashed line alongside a 
fit with the canonical slope of /3 = —2.0. 



/3 = —2.64 ± 0.12 from a much smaller sample of 9 solar 
analogs comparing Lx with P ro t- Our slope is derived 
from a larger sample of stars with a greater range of 
masses and incorporates the convective turnover time to 
compare Rx with Ro. Such a steep slope also supports 
the findings of Fci gelson et al. (20041 ) who found an un- 
expectedly steep decay of X-ray emission as a function of 
age for stars detected in the Chandra Deep Field North, 
which could indicate a steepening of the rotation-activity 
relationship. This slope also produces a much better fit 
to the Sun's Rx, over-predicting it's luminosity by only 
0.2 dex. 

Our fitted slope is inconsistent with the canonical value 
of j3 = —2 to a significance of — 5cr. One of the implica- 
tions of such a steep slope is that these observations are 
inconsistent with the model of a distributed dynamo op- 
erating throughout th e convection zone . the efficiency of 
which scales as Ro~ 2 (jNoves et al. 19841 ). Current mod- 
els of stellar dynamos (at least for solar-type stars with 
radiative cores) use a two-layer interface dyna mo (e.g. 
iParker 1993 : ICharbo nneau fc MacGregor 19971 ). One of 
the implications of such a dynamo is that the dynamo ef- 
ficiency does not scale simply as Ro~ 2 , but is a more com- 
plex function of a number of parameters, not all of which 
can be currently measured or d etermined from models. 

Based on a scaling analysis, Montesin os et al. f 20011 ) 
derived the following expression for the interface dynamo 
number, Nd 



Nn = 



1 

Rd 2 



1 ACl TrrL 



v n d 2 



(2) 



where v is the ratio of the turbulent to magnetic dif- 
fusivities, f2 is the angular velocity at the base of the 
convection zone, AJ7 is the change in angular velocity 
through the tachocline, r cz is the radius of the base of 
the convective zone, d is the characteristics length scale of 
convection at r cz and L is the characteristic length scale 
of differential rotation in the tachocline. We can inter- 
pret the deviation from the proportionality Nr> cx Ro~ 2 



indicated by the unbiased sample of X-ray luminosity ra- 
tios illustrated in Figure [3] in te rms of the additional ex - 
pression in Nd in equations [5] iMontesinos et aL~7 2001) 
argue that the term in brackets can be further simplified 
by appealing to the similar scalings of L and d. Both 
are proportional to the pressure scale height, H p at the 
convection zone base, such that 
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Substituting into Equation [2] we are left with a color- 
dependent term, r cz /H pi and the differential rotation in 
the tachocline expressed as a fraction of the angular ro- 
tation velocity, Afi/fi. The unbiased sample used here 
covers only a small range in effective temperature, there- 
fore for unsaturated magnetic activity it is reasonable to 
assume that r cz /H p does not vary and is also indepen- 
dent of rotation velocity, although these factors might 
be responsible for some of the scatter about the mean 
relation in Figure [3J Recalling that the Rossby number 
is simply Ro — P ro t/T c — 2tt/D,t c , and that t c is approx- 
imately constant for our unbiased sample, the slope of 
the relation in Figure [3] implies that, for solar type stars, 
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i.e. the fractional differential rotation in the unsatu- 
rated regime scales with the angular rotation rate ap- 
proximately according to f2 2 / 3 . The implication of this 
is that the differential rotation of solar-type stars grad- 
ually declines as they spin down. 

3.2. The saturated regime 

Coronal X-ray saturation is known to occur for all 
fast-rotating solar- and late-type stars. Its onset oc- 
curs at larger rotati on periods for stars of lower masses 
(|Stauffer et al. 1994]) . the saturation threshold scaling 
with the convective turnover time, P sat cx r, such that 
it occurs at a const ant value of the Rossby number 
(|Pizzolato et al. 20031 ). The sample used here confirms 
this, showing little variation in Ro sat when the fits de- 
scribed above are performed in different color bins, the 
best fitting value for the entire sample being Ro sa t — 
0.13 ± 0.02. However, if the saturated and unsaturated 
regimes are fit separately (which results in a steeper slope 
to the unsaturated regime) , the interception of these fits 
is shifted to Ro sat = 0.16 ± 0.02. 

That the criteria for saturation is a parameter re- 
lated to the efficiency of the stellar dynamo has 
been argued by some to suggest that coronal satu- 
ration is caused by saturation of the dynamo itself 
(e.g. iGilman 19831 : IVilhu fc Walter 1987| ). a theory sup- 
ported by recent observations of magnetic flux satu- 
ration in rapi dly rotating M dw arfs at Rossby num- 
bers of -0.1 © einers et al. 20091) . However the sat- 
uration thresh olds of other activity indicators are not 
so coherent. Marsdcn et al. (20 09) find saturation of 
the chromospheric Ca n emission lines at Ro — 0.08, 



while iMamaiek fc Hillenbrand (2008L their Figure 15) 
do not find a clear agreement between coronal and 
chromospheric saturation o f field stars. Furthermore 
ICardini fc Cassatella (20071 ) do not find a clear satura- 
tion threshold from Mg n chromospheric emission line 
measurements, and even argue that rotation period is an 
equally good diagnostic for the Mg II line strength. 

Other explanations for X-ray saturation have included 
a saturation o f the filling factor of active regions on the 
star's surface ; Vilhu 19841 ). motivated by a strong corre- 
lation between saturated Lx and stellar radius, rather 
than between Lx and surface temperature. However, 
the detection of rotational modula tion in some satu- 
rated stars (e.g. iGuedel et al. 199 5). and the need for 
enhanced plasma densities to explain t he observed levels 
(|Wood et al. 1994 iDrake et al. 20001 iTesta et al. 2005h 
argues against this theory. Uardine fe Unruh~ 1999) 
suggested that X-ray saturation could alternatively be 
caused by centrifugal stripping of the corona at very fast 
rotation rates. However it will be shown in Section|4]that 
the mass-dependent functional forms of the observed sat- 
uration limit, Ro sa t{M), and the fractional corotation ra- 
dius, r c j i?*(M), are significantly different, which argues 
against such a theory to the first order. 

For our sample we find a mean saturation level of 
Rx = —3.13 ±0.22 (la), almost independent of spectral 
type. This value is in good agreement with the two-part 
power law fit, though it i s less well constrained . As found 
by previous authors (e.g. iPizzolato et al. 20031 ) we find a 
significantly lower mean saturation level for the highest- 
mass stars in our sample (F-type stars), though we will 
show in Section [4] that these stars are not in fact satu- 
rated, but are super-saturated and we therefore conclude 
that the X-ray saturation level is independent of spectral 
type. 

3.2.1. The saturation threshold and the angular momentum 
evolution of stars 

To understand the physical origin of coronal satu- 
ration, the threshold for saturation can be converted 
from one in terms of the Rossby number (though ef- 
fectively the rotation period and the spectral type) to 
age using know ledge of the rotational evolution of stars. 
iBarnes (2003b| ) outlined an empirical formulation of the 
rotation period evolution of solar and late-type stars 
based on open cluster and field star data, character- 
izing two sequences in the P ro t — (B — V) plane: the 
convective (C) sequence of young, fast rotators, and the 
interface (I) sequence of slow rotators. This empirical 
characterization allows ages to be derived from rotation 
periods for stars on the I sequence, a method known as 
gyrochronology. This technique results in a precision of 
±0.05 dex in log r/yr ( excluding absolute uncertaintie s 
in the cluster age scale, IMamaiek fc Hillenbrand "2 008). 
Stars on the C sequence are harder to date accurately 
due to star-to-star variations in the ZAMS arrival time 
and effects such as disk-locking in the pre-main sequence 
p hase. Here the recent gyrochron ology parameterization 
of lMamajek fc Hillenbrand (20081 ) has been used to con- 
vert the empirical threshold for saturation, Ro > 0.13, 
into a threshold in mass-age space, as shown in Fig- 
ure |U It is worth noting that, at the ages considered 
here, d ifferences betwe e n the gyrochronological parame - 
ters of IBarnes (2003bD . IMamaiek fc Hillenbrand (20081) . 
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Fig. 4. — The saturation threshold Ro < 0.13 in mass- 
age space (solid line), compared to the age at which stars 
transition from the rotational C sequence to the I sequence 
(dashed line). All ages calculate d using the rotational evolu- 
tion formulated b y Barnes (20033) an d the parameterization of 
Mamaick & Hillenbrand (2008). 



and lMeibom et al. (20091) are too small to be discernible 
in this figure. 

Also shown in Figure [4] is the mass-dependent age 
at which a star transitions from the C sequence to 
the I sequence. The form of these two, independently- 
determined, empirical transitions is remarkably similar. 
This implies that stars with saturated X-ray emission are 
found on the C sequence of fast rotators, while those with 
non-saturated emission can be found on the I sequence 
of slow rotators. A similar result to t his was already 
hinted at from a much smaller sample by IBarnes (2003al ) 
who associated the presence of stars on the rotational 
I sequence with unsaturated X-ray emitters, but sug- 
gested that stars on the C sequence were coronally super- 
saturated and stars in the gap between the C and I se- 
quences showed saturated X-ray emission. The associ- 
ation of super-saturated coronal emitters with stars on 
the C sequence can be argued against simply based on 
the very small number of super-saturated X-ray emitters 
(see Section ET3]) compared to the much larger number of 
fast rotators found on the C sequence of young clusters. 

The correlation between the transition from saturated 
to unsaturated coronal emission and from the rotational 
C sequence to the I sequence is not surprising when 
one considers that both angular momentum loss and 
coronal X-ray emission are products of the stellar mag- 
netic dynamo and known to scale with th e convective 
turnover time (e .g. iMontesinos et al. 20011 ) . This scal- 
ing was noted bv lBarnes fc Kim (20l6T who showed that 
the functional dependence of the angular momentum 
loss rate was related to the convective turnover time, 
and henceforth derived a new model for the rotational 
evolu tion of cool stars in terms of the Rossby num- 
ber ([Barnes 2010h . In fact the balance between coro- 
nal losses and angular momentum loss is believed to 
be due to the configuration of magnetic field lines at 
the stellar surface, with the fraction of open magnetic 
field lines determining the spin-down rate, while closed 
magnetic field lines dictat e the X-ray luminosity (e.g. 
iHolzwarth fc Jardine 20071) . 
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The similarity between two independently determined 
empirical transition criteria is unlikely to be a coinci- 
dence. It therefore suggests that the changes that occur 
within a star as it transitions from the C sequence to 
the I sequence are also responsible for the star leaving 
the saturated regime of X-ray emission. One could also 
consider the opposite causality, that the mechanisms re- 
sponsible for the changes in X-ray emission also cause 
changes in the angular momentum loss rate. However, 
the evolution of angular momentum loss is a more com- 
plex subject that must explain other observable features 
(e.g. the sequence of young ultrafast rotators in clus- 
ters) and it therefore seems unlikely that the mechanism 
responsible for X-ray saturation is responsible for these 
effects as well. 

The hypothesis that we are therefore left with is that 
the processes governing stellar angular momentum loss 
are also responsible for X-ray saturation. iBarnes (2 003b) 
argued that the C and I sequences are due to the coupling 
of the stellar wind to, respectively, just the convective 
zone (which is decoupled from the radiative zone), and 
to the entire star. The transition between these zones 
(across the 'gap' in the color-rotation period diagram) 
is therefore associated with the couplin g of the radia - 
tive and convective zones to each other. IBarnes (2 003b) 
suggests that stars on the convective sequence generate a 
convective or turbulent dynamo (e.g. lDurnev et al. 19931) 
and that the shear between the fast spinning radiative in- 
terior and the convective envelope eventually generates 
an interface dy namo that results in the transition onto 
the I sequence. IBarnes fc Kim (20101 ) have since argued 
that the form of the rotational isochrones does not agree 
with the form expected for the relevant moments of in- 
ertia outlined in this model. However, the association 
of the two regimes with two different dynamos is worth 
exploring further. 

Here we suggest that coronal saturation does not re- 
quire a separate physical mechanism, but is a manifes- 
tation of the different dynamos that are present in stars 
on the convective and interface sequences. In this pic- 
ture, slowly rotating stars have a cyclic, interface-type 
dynamo with Skumanich-style spin-down, while faster ro- 
tators have a turbulent dynamo. The efficiencies of the 
two different dynamos could have significantly different 
dependencies, in terms of stellar parameters, and there- 
fore have different forms in the Rx~Ro diagram. X-ray 
saturation would therefore not be an actual saturation in 
any sense, but a completely different magnetic dynamo 
configuration with completely different dependencies. 

If this scen ario is correct it partl y answers the ques- 
tion raised bv lPizzolato et al. ( 20031) as to why the bolo- 
metric luminosity is the only parameter necessary to de- 
termine the X-ray luminosity in the saturated regime, 
despite the Rossby number being considered the funda- 
mental parameter of the efficiency of the stellar dynamo. 
The answer is that there are two different dynamo con- 
figurations, a convective dynamo parameterized by the 
bolometric luminosity, an interface dynamo parameter- 
ized by the Rossby number. This does however raise the 
question of why the transition between the two regimes 
is apparently so smooth in the Lx/L(>oi versus Ro dia- 
gram. In other words, why is the luminosity ratio for 
unsaturated stars close to the saturation threshold equal 
to ~ 10~ 3 ? The answer to this question will require a 



much better understanding of the two dynamo configura- 
tions present in the saturated and unsaturated regimes. 

3.3. Supersaturation 

At very high rotational velocities the fractional X- 
ray luminosity h as been observed to decrease below the 
saturation level (jRandich et al. 19961 ) . an effect dubbed 
"supersaturation". This effect has only been persua- 
sively observed in the coronal emission from young (< 
100 Myrs) G and K dwarfs and is primarily based 
on projected rotational velocities, not rotation peri- 
ods, in the y oung clusters IC 2391, IC 2602, and 
a Persei (e.g. Randic h et al. 19 96: Pro sser et al. 19961 : 
Stauffer et al. 19971: IStepieh et al. 20011 ) from which 
Stauffer et al. (19971 ) found a decline in X-ray satura- 
tion for stars with Ro < 0.01. Despite this, evidence 
for supersaturation in l ower-mass stars is s carce; both 
I James et al. (20001 ) and Ueffries et al. (20111 ) could find 
no evidence in samples of ultra-fast rotating M dwarfs. 

Like coronal saturation, the phenomenon of super- 
saturation and its physical origin is still heavily debated. 
Figure[5]shows Lx / L^oi for all stars with saturated X-ray 
emission as a function of various quantities suggested as 
parameters influencing the supersaturation effect. The 
first two panels of Figure O show the effect of rotation 
period and Rossby number on the X-ray luminosity ra- 
tio, neither of which show evidence for supersaturation, 
or a notable decline in the luminosity ratio for either the 
fastest rotators or the sources with the smallest Rossby 
n umbers. 

Uames et al. (2000) developed the theory of 

Jardi ne fc Unruh (19991 ) that the X-ray emitting 
volume of rapid rotators is reduced via centrifugal 
stripping. This theory was initially put forward to 
explain X-ray saturation, with the increase in coronal 
density balancin g the decrease in v olume and leading to 
saturation, but Uames et al. (2000l ) and Uardine ( 2004) 
used it to explain super-saturation for stars rotating fast 
enough that the Keplerian corotation radius becomes 
very close to the stellar surface. As a function of 
the stellar radius , the corotation radius is given by 
Uames et al. (20001 ) 

RKepler/R* = (^^) V3 « P 2/3 (6) 

which, using stell ar properties taken from the models of 
ISiess et al. (200(3 ). allows the fractional X-ray luminosity 
to be calculated as a function of the relative corotation 
radius, as shown in Figure [S] In this sample there ap- 
pears to be some evidence for a decrease in the typical 
fractional X-ray luminosity of stars with the lowest coro- 
tation radii. Based on the mean Lx/Lboi a decrease is 
se en for RKevier / 1 R* 5 , 3, in agreement with t hat found 
bv l James et al. (20001) and lJeffries et al. (201 1[ ). This ef- 
fect is still evident if one considers only K stars, and to 
a lesser extent G stars, but not M stars, as found by 
previous authors. As parameterized in Figure [5] the bin 
containing the fastest rotators hints at a moderate de- 
cline to Lx/Li, i = —3.41 ± 0.26, deviating la below 
the mean saturation value. A two-sample Kolmogorov- 
Smirnov (KS) test comparing the X-ray luminosity ratio 
of stars above and below the proposed corotation super- 
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Fig. 5. — X-ray to bolometric luminosity ratios for all stars with saturated X-ray emission (Ro < Ro sa t) as a function of the rotation 
period (upper left), the Rossby number (upper right), the Keplerian co-rotation radius (lower left), and the excess polar updraft acceleration 
(lower right). Stars are colored as per their spectral type (F-type stars are blue, G-type stars are green, K-type stars are black, and M-type 
stars are red). Also shown are 25-star average luminosity ratios, with standard error and standard deviation error bars in full and dotted 
lines respectively. 



saturation threshold of RKepler / ' R* ~ 3 gives a probabil- 
ity of < 10~ 3 that the two samples are drawn from the 
s ame population. 

iStepieri et al. (20011 ) offered an alternative explanation 
that the reduction in X-ray luminosity is caused by a de- 
crease in the filling factor of coronal active regions on the 
stellar surface. They propose that a poleward migration 
of active regions would be induced by nonuniform heat- 
ing of the conve ctive envelope in a ccordance with the von 
Zeipel theorem. Ivon Zeipel (1924) argued that the effect 
of rotation on the internal radiative transfer within a star 
causes a higher fraction of flux at the poles than at the 
equator, an effect known as gravity-darkening. For suf- 
ficiently fast rotating stars this flux imbalance can also 
induce str ong convective updr afts in the outer convective 
envelope (jStepien et al. 19981) , capable of sweeping mag- 
netic flux tubes in the interior to the poles before they 
rise to the surface. 

According to the von Zeipel theorem the emergent sur- 
face flux for a rotating solid body is proportional to the 
local effective gravity, g e ff, given by 



9eff(u) 



GM 



u 2 Rsin 2 6 



(7) 



where w is the angular rotation rate, R is the radius, 
and 9 is the colatitude (differences in t he theorem fo r 
a differentially rota ting star are minor, iMaeder 1 999). 
IStepieh et al. (20011 ) have argued that the resulting 
change in X-ray luminosity is proportional to the relative 
strength of this "polar updraft" effect at the boundary 
between the radiative interior and the convective enve- 
lope, which we quantify as 



Gx = 



gef/M -g(0) 
5(0) 



w 2 R c sin 2 
GM/R 2 



(8) 



where R c is the radius of the radiative-convective bound- 
ary. The efficiency of this mechanism therefore scales as 
the rotation rate to the second power and the radius of 
the radiative-convective boundary to the third power. 
A dopting the radiu s of this boundary from the models 
of iSiess et al. (20"00l as a function of both spectral type 
and age), Figure [5] shows the X-ray luminosity ratio as a 
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function of the absolute value of this parameter, with a 
clear decline for Gx > 0.01. The bin containing the stars 
with the strongest polar updraft shows a deviation from 
the mean level with a value of LxjLboi = —3.45 ± 0.16, 
approximately 2a below the mean value for the satu- 
rated stars 6 . A KS test comparing stars above and below 
the proposed polar updraft supersaturation level gives a 
probability of < 10~ 3 that the two samples are drawn 
from the same population. 

Based on the sample analyzed here, while no single 
theory provides a conclusive description of the super- 
saturation phenomenon, we can show that neither the 
rotation period nor the Rossby number are an accurate 
diagnostic of the effect, while either the Keplerian coro- 
tation radius or the excess polar updraft offer a clearer 
trend. Of these two, Stepien's theory of poleward migra- 
tion of active coronal regions has the stronger correlation 
with the observed decrease in the X-ray luminosity ratio 
and a more significant deviation from the mean X-ray 
luminosity ratio saturation level. 

The fundamental difference between these two theories 
for supersaturation concerns the underlying structure of 
the coronal loops that give rise to X-ray emission. Coro- 
nal stripping will occur if the loops are large, low-density, 
and unstable to the centrifugal force. In Figure [5] a de- 
crease in the luminosity ratio begins at Rxepier / R* ~ 3, 
which would suggest that coronal loops are as large as 
two stellar radii in height. Alternatively, coronal loops 
may be compact, high density, and influenced by the dy- 
namics of surface flows that induce a polar updraft of 
coronal regions towards the stellar poles. The difference 
between the two supersaturation theories is therefore the 
size and density of the coronal plasma. High-resolution 
X-ray spectroscopy (necessary to determine such quanti- 
ties) of nearby, supersaturated stars are very rare. The 
only such observation to date is that of the coronally 
supersaturated (log Lx/Lboi = —3.6, Ro ~ 0.03) star 
VW Cephei by iHuenemoerder et al. f 20061 ) who found a 
coronal density of log N e = 10.5 — 11.25. For VW Cephei 
this is equivalent to a coronal loop height of 0.06-0.2 stel- 
lar radii, significantly below that at which supersatura- 
tion would be expected to occur (Figure [5| or the coro- 
tational radius of VW Cephei itself, Rxepier/R* — 1-75. 
These observations therefore suggest that coronal loops 
are compact and high-density and sufficiently far below 
the corotation radius of the star that coronal stripping 
is unlikely to be responsible for supersaturation. Simi- 
lar coronal densities have been o bserved for a sam ple of 
saturated and unsaturated stars (|Testa et al. 20041 ). 

While it is not clear that the heating mechanism 
for the chromosphere is the same as those for the 
corona, their emission scales reasonably well and there- 
fore chromospheric emission can offer another perspec- 
tive on stellar activity. If coronal supersaturation 
were caused by the migration of active regions towards 
the poles one might expect to s ee supersaturation in 
chr omospheric emission as w ell. iMarsden et al. (20091 ) 
and iJackson fc Jeffries (20101 ) have searched for chromo- 
spheric supersaturation in a sample of F, G, K, and M- 

6 The variance in the luminosity ratio for the supersaturated 
stars under the parameterization of polar updraft is smaller than 
that under coronal stripping, as is clear from Figure[5] This results 
in a more significant deviation of the supersaturated stars from the 
mean value under this parameterization. 



type stars but could find no evidence for a reduction 
in chromospheric activity indicators from a small sam- 
ple of ultra-fast rotators. Since chromospherically active 
regions with their lower scale height should not be vul- 
nerable to centrifugal stripping in the same way coronal 
active regions are, the lack of evidence for chromospheric 
supersaturation, if confirmed, would argue for coronal 
stripping as the likely mechanism for supersaturation. 

The mechanism behind supersaturation remains a mys- 
tery despite this work. Coronal and chromospheric ob- 
servations of even faster rotating stars that might exhibit 
stronger levels of supersaturation are necessary to resolve 
this. Such observations could provide a more significant 
proof of the lack of chromospheric supersaturation, or ev- 
idence for chromospheric saturation at a faster rotation 
rate. Increasing the sample of supersaturated stars would 
also allow the correlations studied here to be probed in 
greater detail. 

4. CORONAL X-RAY EMISSION OVER STELLAR 
LIFETIMES 

In this section an attempt is made to combine the re- 
lationships and thresholds determined in this work to 
broadly sketch out the changes in coronal X-ray emis- 
sion that occur over stellar lifetimes, as a function of the 
rotation period. Figure |6] shows the mass-rotation plane 
for all stars in our sample. Also shown are the borders of 
the three regimes of unsaturated, saturated, and super- 
saturated coronal X-ray emission. The super-saturated 
regime is defined both using the theory of coronal strip- 
ping (assuming a threshold of r c /i?* ~ 3) and the theory 
for polar updraft (using Gx ~ 0.01). The similarity be- 
tween the theoretical predictions of the supersaturation 
thresholds of these two theories is particularly notable, 
despite their different theoretical formulations. This goes 
a long way towards explaining why both parameteriza- 
tions in Figure [5] showed significant evidence for the su- 
persaturation effect. The forms of both of the supersat- 
uration criteria also illustrates why only G and K dwarfs 
have been observed to have super-saturated X-ray emis- 
sion, since M dwarfs would have to be rotating excep- 
tionally fast to exhibit the phenomenon, possibly at the 
break-up period, and certainly faster than the typically 
observed rotation periods for young, ultrafast rotators of 
-0.2 days. 

Also clear from Figure [6] is that the majority of F- 
type stars (M > 1.15M Q ) will never exhibit satu- 
rated X-ray emission over their lifetimes, passing from 
the super-saturated regime straight to the unsaturated 
regime. This exp lains why previous authors (e.g. 
iPizzolato et al. 20031 ). as well as this study (see Sec- 
tion [5]), find a lower X-ray luminosity ratio for saturated 
F-type stars, (L x / Lboi) sat — —4.3, than for solar and 
lower-mass stars, (Lx / Lboi)sat ~ —3.2. Since they are 
not truly in the saturated regime their emission should 
not be treated as such. 

The theoretical form of constant corotation radius in 
the mass-rotation diagram is also notably different from 
the empirically-derived form of the saturation boundary, 
Ro — 0.13, which suggests that centrifugal stripping is 
unlikely to be responsible for X-ray saturation 7 . 

7 A caveat to this is that X-ray saturation via centrifugal strip- 
ping is suggested to be a balance between the decreasing coro- 
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TABLE 2 

Empirical convective turnover times 



V - K s B — V Mass (M ) Fit with g = -2 Fit with g = -2.70 

range range range log H X sat F aat (days) log r (days) ~ log % sat fsat (days) log r (daysj 
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Empirical convective turnover times fitted from the data divided into V — K s color bins (with equivalent B — V color and main sequence 
mass bins also shown). Fits were performed assuming a constant saturation limit, Ro = 0.13, as determined above, and the slope of the 
unsaturated rotation - activity regime to be either /3 = —2 (the standard value) or /3 = —2.70, as estimated from the fit to to our 
unbiased sample, lcr uncertainties are also shown, estimated using a bootstrapping technique. 
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Fig. 6. — The three regimes of coronal X-ray emission shown in 
Mass-P ro t space, with all sources in this sample shown as black 
dots. The transition from saturated to unsaturated emission is 
defined by Ro > 0.13 (blue dashed line). The transition from sat- 
urated to supersaturated emission is shown for both the theories 
of coronal stripping (r c /ij* ~ 3, red dashed line) and coronal up- 
drafts (Gx ~ 0.01, green dashed line). 



5. EMPIRICAL DETERMINATION OF THE CONVECTIVE 
TURNOVER TIME 

It is clear from the work performed here that a 
fundamental quantity in the rotation-activity rela- 
tion and the evolution of angular momentum loss 
for all cool stars is the convective turnover time, r 
(e.g. Figure [2]). This color- dependent quantity can 
be estimated theore tically (e.g. iKim fc Demarque 19961 : 
Ventura et al. 19981) or determ ined empirically (e.g. 
Noves et al. 1984 IStepien 19941 ). though both methods 
have mainly been applied to sola r-type stars, with very 
little work on very low mass stars. iPizzolato et al. (20031 ) 
were the first to extend the empirical approach into the 
M-dwarf regime, but even their study was limited to only 
21 stars with B — V > 1.55. The sample used here is al- 

nal volume (as the corotation radius decreases) and an increasing 
plasma density (as rotation rate i ncreases) that results in an ap- 
parently uniform saturation level (Jardinc & Unruh 1999), there- 
fore the comparison of these two curves does not represent the 
complete picture. 



most an order of ma gnitude larger than that used by 
IPizzolato et al. (20031 ). which contained 164 stars with 
B -V > 1.55 and 26 stars with B -V > 1.7, and could 
therefore offer a considerably improved empirical deter- 
mination of t\ 

Following IStepien (1991 we derive an empirical esti- 
mate of t by assuming that the Lxj^boi — Re diagram 
can be divided into saturated and non-saturated regimes 
with the formulations described above. To do this we 
have divided our sample into 10 bins in V — K s or B — V 
color, with an approximately equal number of stars in 
each bin (see Table \2§ and then fitted the equations 

L X \ C B -V Prot forP™* > Psat 

^ iteL ^P rot <P sat 

varying the parameters (Lx I Lboi)sat and P sa t to provide 
the best fit. This was performed using j3 values of -2.0 
(the canonical value) and -2.70 (the value determined 
here). From these fits the color-dependent constant, 
Cv-K 3 = Ct 2 , may be determined and, therefore, the 
color-dependence of the convective turnover time, r, by 
setting the scaling constant, C, so that our values of r for 
solar-mass stars match those of iNoves et al. (T984[) for 
the Sun. The results of this fitting process and the em- 
pirical determinations of the convective turnover time are 
listed in Table [2] with la confidence intervals determined 
using a bootstrapping technique. Also shown are the 
main-sequence masses of each bin, converting the V — 
color to mass using the models of lSiess et al. (20"00l) . 

Figure [7] shows a comparison of our values with pre- 
vious empirical estimates from the literature. Our new 
empirical values are in good agreement with previous es- 
timates but have smaller error bars. At the high-mass 
end we do not observe the tu rnover towards lowe r convec- 
tive turnover times found bv lNoves et al. ( 19841 ) . though 
their data only extended to B — V ~ 0.45 and are ex- 
trapolated beyond that. As shown in Section |4j the fact 
that late F-type stars may never pass through the satu- 
rated regime will affect the determination of convective 
turnover times with this method and so also the validity 
of both our results and previous methods that ignored 
the effects of supersaturation. 

At low masses our results show a rise in the convective 
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in log r of 0.031 dex, an improvement of a factor of two 
over the single power-law fit. 

The different appearance between the two parameter- 
izations as a function of color in Figure [7] is almost cer- 
tainly due to the different dependencies of the two col- 
ors on whatever is the relevant physical parameter that 
determines the convective turnover time. In Figure [7] 
we also show the convective turnover time as a function 
of main-sequenc e stellar mass us ing the 1 Gyr stellar 
isochrones from ISiess et al. (20001 ). Under this param- 
eterization the convective turnover time shows a much 
smoother rise towards lower masses with no evidence 
of the discontinuity seen as a function of the two col- 
ors. This could suggest that the stellar mass is the rele- 
vant physical parameter that determines the convective 
turnover time, or at least is more closely related to it 
than the two colors that have complex and non-linear 
dependencies on most stellar parameters. 

An OLS (Y|X) fit gives the linear relation logr = 
2.12 - 0.96 (M/M Q ), with an rms of 0.060 dex, but a 
much better fit can be obtained using a second order 
log-log polynomial 



Mass / M e 

Fig. 7. — The empirically-determined convective turnover time 
as a function of V — K a color (top), B — V color (middle), and main- 
sequence mass (bottom), using the best-fitting model determined 
above (/3 = —2.70) and with lcr error bars for all points. In the 
top panel the two-part power-law fit described in the text is shown 
with a red line. In the middle p anel our values are compared with 
the empirical determi nations of [Noycs ct al. (1984], blue line) and 
Pizzol ato et al. (20031 . red points). In the bottom panel a second 
order polynomial fit as a function of log (M) is also shown. Note 
that the x-axis of the bottom panel has been reversed to mimic the 
color-scale of the other two panels. 



turnover time as found by iPizzolato et al. (20 03) . This 
rise occurs around B — V ~ 1.4, or spectral type M0, 
significantly before the fully-convective boundary that 
occurs at types M3-M4 (B — V ~ 1.5-1.7). Unlike 
previous studies we can resolve this rise and find that it 
begins to plateau again around B — V ~ 1.6, which is 
closer to the point at which stars are believed to become 
entirely convective and where it has been suggested a 
different dynamo mechanism may be at work. However, 
if a different dynamo mechanism is at work, then the 
parameterization used here to calculate the convective 
turnover time may not be accurate and therefore these 
values, and the entire use of the convective turnover time 
for such stars, may be invalid. 

The rise in r seen at B — V ~ 1.4 is less evident when 
one considers the convective turnover time as a function 
of the V — K s color (Figure [7]) , where the relationship 
between the two quantities is almost exponential. An 
OLS (Y|X) fit to log t as a function of V — K s provides a 
simple relation logr = 0.64 + 0.27 (V — K s ) with an rms 
dispersion in r of 0.058 dex. However a discontinuity is 
still visible at V — K s ~ 3.75 (type M0) and therefore a 
two-part power-law is more appropriate. This was fit as 



lo 



0.73 + 0.22 A 
-2.16 + 1.50 A 



forA < 3.5 , , 
0.13 A 2 forA > 3.5 [W) 



where A = V — K s ; the fit is valid over the range 1 . 1 < 
V — K s < 6.6. This two-part fit has an rms dispersion 



logr = 1.16- 1.491og(M/M ) - 0.541og 2 (M/M Q ) (11) 

which has an rms dispersion in log r of 0.028, the lowest 
of any of the fits here, and valid over the range 0.09 < 
M/Mq < 1.36. 

It should be noted that the validity of any of these 
relationships are still highly uncertain for V — K s > 5 
(B — V > 1.6, M < 0.15M©) due to the paucity of data 
for low-mass stars. Despite this, these are certainly the 
most accurate and comprehensive empirical estimates of 
the convective turnover time for solar- and late-type stars 
to date. Observations of slowly rotating, X-ray unsatu- 
rated M dwarfs will be necessary to improve these but 
are, unfortunately, very difficult to obtain (but see e.g. 
iRobrade fe Schmitt 20091 ). 

6. CONCLUSIONS 

The largest sample to date of solar and late-type stars 
with well measured X-ray luminosities and photometric 
rotation periods has been gathered from the literature. 
This sample of 824 stars is over three times larger than 
previous samples and is used here to study the relation 
between rotation (in the form of the Rossby number, 
Ro = P rot I t) and stellar activity (in the form of the X- 
ray luminosity ratio, Rx = Lx/Lb i). The Rx — Ro 
plane is divided into the unsaturated, saturated, and 
super-saturated regimes of coronal X-ray emission, and 
a number of fits are performed on these regimes. 

In the unsaturated regime, where the relationship be- 
tween Ro and Rx probes the efficiency of the stellar 
dynamo, we use a smaller unbiased sample of 36 stars 
to fit a slope of f3 = —2.70 ± 0.13, significantly steeper 
than the previously assumed value of j3 = —2. Such 
a steep slope is inconsistent with a distributed dynamo 
model and lends we ight to the theory of an interface dy- 
namo (|Parker 1993f ) and an additional term in the dy- 
namo number equation . Using the scaling analysis of 
iMontesinos et al. (20011 ) for the interface dynamo num- 
ber, we show that this result implies Af2/f2 oc fi - 63 for 
solar-type stars, i.e., that the fractional differential rota- 
tion in the unsaturated regime scales with the angular 
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rotation rate to the power of two thirds. The steeper 
relationship between rotation rate and activity level also 
provides a better fit to the Sun's activity level than the 
canonical (3 = —2 slope that over-fitted the Sun's X-ray 
luminosity by an order of magnitude. 

In the saturated regime we find a mean saturation level 
of Rx = —3.13 ± 0.22, which we argue is independent of 
spectral type once the super-saturated F-type stars are 
removed. We explore evidence for the super-saturation 
phenomenon by parameterizing the activity level as a 
function of a number of parameters suspected of causing 
supersaturation. We find that the parameters of rotation 
period and Rossby number do not show a significant cor- 
relation with activity level, while the parameters of the 
corotation radius and the excess polar updraft both show 
a significant deviation from the mean saturation level. 
Two sample KS tests comparing saturated and super- 
saturated stars (using both parameters to separate the 
two populations) find a negligible probability that they 
are drawn from the same population. The the ory of po- 
lar updraft suggested by IStepieri et al. C2001D is found 
to provide the most significant deviation from the mean 
saturation level, and combined with other observations 
we argue that this is the more likely explanation for the 
supersaturation phenomenon. 

The color-dependent form of both supersaturation the- 
ories are shown to be significantly different from the 
empirically-determined saturation criteria, such that it 
is unlikely that either are responsible for the satura- 
tion of coronal emission, which we find to occur at a 
Rossby number of Ro = 0.13. Instead, based on the co- 
incidence between the empirically determined saturation 
threshold and the time at which stars transition from the 
rotational cony ective sequence to the interface sequence 
(|Barnes 2003bD , we suggest that the phenomenon of dy- 
namo saturation is actually caused by a change in the dy- 



namo configuration of the star. This theory precludes the 
need for an actual saturation effect and argues that the 
two main regimes of coronal emission should be treated 
as separate probes of different dynamo configurations. 

Finally this sample is used to determine a new estimate 
of the empirical convective turnover time as a function of 
the B — V color, the more preferable V — K s color, and 
the main-sequence stellar mass. When parameterized as 
a function of either photometric B — V color, the con- 
vective turnover time shows an increase around spectral- 
type MO, notably before the fully convective boundary 
(~M3-M4). A two-part power- law fit to the convective 
turnover time as a function of V — K s color provides the 
simplest method for estimating this parameter with read- 
ily available photometry. When parameterized as a func- 
tion of stellar mass, the convective turnover time shows 
a much smoother rise towards lower masses with no evi- 
dence for a discontinuity, possibly hinting at the relevant 
physical parameter. The convective turnover time can 
be well fit as a second order log-log polynomial. 
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APPENDIX 

APPENDIX A : OBSERVATIONS OF FEPS STARS 

We acquired time-series photometric observations of 100 stars drawn f rom the Spitzer Legacy Science Program For- 
mation and Evolution of Planetary Systems (FEPS, iMever et al. 2 006). Nightly photometric observations of these 
stars were acquired with several of the Automatic Photometric Telescopes (APTs) operated by Tennessee State Uni- 
versity (TSU) and located at Fairborn Observatory in southern Arizona. The telescopes included the T2 0.25 m, the 
T3 0.40 m, the T4 0.75 m, and the T10, Til, and T12 0.80 m APTs. Each FEPS star was observed differentially with 
respect to two or three comparison stars; all observations were made through standard filters and transformed to the 
Johnson BVRI or to the Stromgren by photometric system. Further information about the design and operation of 
th e APTs as wel l as the o bser ving and data reduc tion techniques and the resulting photometric precision can be found 
in lHenrv Vl99^allbiri999l) and lEaton et al. (20031) . 

Most of the FEPS stars were observed once each clear night throughout one of their individual observing seasons 
between 2002 September and 2004 July. A few of the stars have been observed for a decade or more bec ause they are 
part of an ongoi ng, long-term program to monitor subtle brightness changes in solar-type stars (see, e.g.. iHenrv 19991 
Hall et al. 20091) . All data sets were subjected to periodogram analysis based on least-squares fitting of sine curves 
over a range of trial periods. Definitive photometric periods were determined for 72 of the 100 FEPS targets, of 
which 65 also have well-constrained X-ray luminosities and appear in our catalog. We take each photometric period 
to be the stellar r otation period made visible by rotational modulation of the visibility of photospheric starspots. 
IHenrv et al. (19 95) provide several dozen examples of low-amplitude light curves of active cool stars discovered with 
these same techniques. 

APPENDIX B : THE ACTIVITY-ROTATION CATALOG, NAMING NOMENCLATURE AND REFERENCES 

The final catalog of 824 sources with X-ray luminosities and rotation periods is available online. The catalo g includes 
a number of stellar parameters derived from the available photometry and the models of ISiess et al. (20"OOT ). For the 
cluster stars these are dependent on the properties of the cluster itself (see Table [I| . To allow the catalog to be easily 
cross-referenced with other works, source names are provided for sources, where available, from the Henry Draper 
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TABLE 3 

List of fields in the Activity-Rotation catalog 



Label 


Units 


Explanation 


Name 




Default source name 


HD number 




Henry Draper catalog number 


HIP number 




Hipparcos catalog number 


GJ number 




Catalog of Nearby Stars number 


GCVS name 




General Catalog of Variable Stars name 


2MASS name 




2MASS source name 


RA 


deg 


Right ascension (J2000) 


Dec 


deg 


Declination (J2000) 


Cluster 




Name of cluster source belongs to (see Table [TJ 


Distance 


pc 


Distance 


V 


mag 


Cousins V-band magnitude 


V - K s 


mag 


Dereddened (V — K s ) color 


SpT 




Spectral type (if used to infer a dereddened V — K s color) 


r_SpT 




Reference for spectral type (see list below) 


log L x 


erg s — 1 


Logarithm of the X-ray luminosity 


r.L x 




Reference for the X-ray data (see list below) 


Prot 


days 


Rotation period 


IC-Prot 




Reference for the rotation period (see list below) 




M Q 


Stellar mass from Siess models 


Ft* 


Rq 


Stellar radius from Siess models 




Lq 


Logarithm of stellar luminosity from Siess models 


T eff 


K 


Stellar effective temperature 


dcz/R 




Height of radiative-convective boundary from Siess models (as a function of the stellar radius) 


log L x /L bol 




Logarithm of X-ray to bolometric luminosity 



Refere nces for spectral ty pes: B46 = IBinnendiik (194j j), B84 =IBreger (19841), B85 = Bidclman (198511 , B 95 =IBaIiunas et al. (19951 ). 
B97 = IBouvier et al. (1997ft, C18 = ICannon & Pickering (191811. Dll = IDelorme et al. (20111), G01 = IGrav et al. (20011), G03 = 
IGrav et al. (20031), G06 = IGrav et al. (2009), H56 = IHeckmann et al (1959) , H62 = IHerbig (196%, H99~= IHouk fe Swift (19991) , 
H02 = IHenrv et al. (200% J74 = Uov fe Abt (19741). K67 = IKraft (19671) . K89 = IKeenan fe McNeil (19891) . L82 = ILu (19821 ). 
M65 = IMorgan fe Hiltner (19651) , M81 = IMundt fc Bastian (19811 ). R41 = IRamberg ( 19411) . R6 6 = IRebeirot (1969 ), S51 = 
I Schwassmann fc van R hiin (1951). S86 = IStephenson (19861 ). T06 = ITorres et al. (2009 ), U72 = lUpgren et al. (19721) . W48 = 

I Wilson (19481) . W63 = I Wilson (1963D . 

Rerences for X-r ay lum inosities: B03 = IBriggs fc Pve (20031). F03 = IFranciosini et al. (20011), J0 3 = Uudge et al. (20031). Jll = 
Ueffries et al. (20111). M90 = lM"icela et al. (19901). M96 = IMicela et al. (199ef). M99 =IMicela et al. (19991). P03 = IPizzolato et al. (20031) . 
P06 = IPillitteri et al. (2009 ). R95 = IRandich fc Schmitt (19951 ). S94 = IStauffer et al. (19941 ). S04 = ISchmitt fc Liefke (20041 ). V99 = 

IVoges et al. ( 19991 . ROSAT all-sky bright sourc e catalog). 

R eferences for rotati on perio ds: B96 = IBaliunas et al. (1999 ). B97 = IBouvier et al. (19971) . C99 = ICutispoto et al. (1999T) . D96 
= IDonahue et al (1991) . Dll = IDelorme et al. (20111) FEPS = FEPS Program (see Appen dix), GE 06 = IGe et al. (2006T) ~H95 = 
I Henrv et al. (19951). H OP = IHenrv et al. (2000) . H06 = IHebrard fc Lecavelier Pes Etanes (20091. H10 = IHartman et al. (20101). Hll = 
lHartman et al. (20111). 107 =IIrwinet al. (20071). 108 =IIrwin et al. (20081) , K02 = IKoen fc Ever (20021) . K07 =IKiraga fc Stepien (20071 ). 
M1 = IMessina et al. (20101 N07 = INorton et al. (20071) . P03 = IPizzolato et al. (20031) . P04 = IPaulson et al. (200*41 ). P05 
= IPoimanski et al. (20051 ). O 01 = IQueloz et al. (20011). S00 = IStrassmeier et al. (20001) . Sll = IScholz et al. (20111) . W98 = 
IWichmann et al. (19981) 7X06 = IXing et al. (2009 ). X07 = King et al. (2007fT . 

catalog (HP. ICannon fc Pickering 19181 186 sources), the Hipparcos catalog (HIP. IPerrvman et al. 19971 22 sources), 
the Catalog of Ne arby Stars (GJ, iGliese fc Jahreifi 1991L 37 sources), the General Catalog of Variable Stars (GCVS, 
Version 2011 Jan. ISamus et al. 2009L 293 sources), and the 2MASS catalog (435 sources). The default name is taken 
from the X-ray or rotation period paper in the case of cluster stars, and for fields stars is chosen according to the 
catalog order listed above. 202 sources remain without a name. 

The final catalog will be made available online. For convenience we include here only a list of the different fields in 
the catalog (Tabled]) and an explanation of their meanings. References for all the X-ray luminosities, rotation periods, 
and spectral types used in this work are also listed as part of Table [3] 
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